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Abstract 

Quark Gluon Plasma transport coefficients for heavy quarks and QQ pairs are computed 
through an extension of the results obtained for a hot QED plasma by describing the 
heavy-quark propagation in the eikonal approximation and by weighting the gauge field 
configurations with the Hard Thermal Loop effective action, ft is shown that such a model 
allows to correctly reproduce, at leading logarithmic accuracy, the results obtained by 
other independent approaches. The results are then inserted into a relativistic Langevin 
equation allowing to follow the evolution of the heavy-quark momentum spectra. Our 
numerical findings are also compared with the ones obtained in a strongly-coupled sce- 
nario, namely with the transport coefficients predicted (though with some limitations 
and ambiguities) by the AdS/CFT correspondence. 
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1. Introduction 

Ultra-relativistic heavy-ion experiments at the Relativistic Heavy Ion Collider (RHIC) 
and (in the near future) at the Large Hadron Collider (LHC) aim at reproducing, in a 
small region of space-time, the conditions of the primordial universe, with such a large 
energy-density to allow the onset of the deconfined phase of QCD. Through high-energy 
nucleus-nucleus collisions one should be able to explore the QCD phase-diagram in the 
region of high temperature and (almost) vanishing baryon density. 

However the fireball of quarks and gluons possibly produced in such collisions expands 
and cools, so that the only strongly-interacting particles reaching the detectors are again 
color-singlet hadrons. Hence the need to extract information on what occurred before 
hadronization: whether a thermalized Quark Gluon Plasma (QGP) was produced and, 
if so, which are its properties. 

Heavy quarks (charm and bottom) and quarkonia are well suited for this purpose. 
Due to their large mass they are produced in the very early stages of the collision and, 
before giving rise to experimentally detectable signals from their decays into hadrons 
and/or leptons, they have to cross the hot (possibly) deconfined region. By comparing 
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the data collected in nucleus-nucleus collisions with the ones coming from benchmark 
experiments, like proton-proton or proton-nucleus collisions, in which one does not expect 
QGP formation, one eventually can prove the formation of a totally new kind of medium. 
In particular, a lot of interest on quarkonia was triggered by the seminal paper by Matsui 
and Satz [T| proposing the anomalous suppression of the J/'J' in nucleus-nucleus collisions 
as an unambiguous signature of deconfinement. 

For what concerns heavy quarks — giving rise to D and B mesons at hadronization 
— they can provide useful insight into the transport properties of the medium they cross. 
Due to their large mass they require a much longer time (by a factor M /T) than the light 
particles to thermalize; naively one expects that in an expanding medium with a quite 
short life-time (of order of 5 — 10 fm/c) they can hardly approach thermal equilibrium. 
If, on the other hand, one found them to follow the flow of the fluid, this would be a 
signal of large friction and momentum-diffusion coefficients, whose origin would need an 
explanation. 

RHIC data led to a description of the matter produced in the Au-Au collisions as a 
strongly-coupled QGP (sQGP). In particular the radial and elliptic (in semi-central 
collisions) ffow [3, (B| observed in the pT-spectra of detected hadrons suggested a picture 
of the QGP in the regime of temperatures accessible at RHIC more similar to a ffuid 
(hence with very small mean-free-paths) — with a collective expansion driven by pressure 
gradients — than to a weakly-interacting gas of quarks and gluons. Furthermore, the 
strong suppression (by a factor of 5) of high-py hadrons [6, 7] with respects to elementary 
proton-proton collisions suggests that the QGP produced at RHIC is a very opaque 
medium. 

Heavy quarks at RHIC are studied through the electrons arising from the semi- 
leptonic decays of D and B mesons. This makes it hard to draw definite conclusions on 
their flow and energy loss, due in particular to the difficulty of disentangling the charm 
and bottom contributions. However, recent measurements performed by the PHENIX Q 
and STAR Q collaborations seem to support a scenario with the charm quark following 
the flow of the medium and experiencing a strong energy loss, comparable to the one 
seen in light hadron spectra. Different explanations where advocated to account for such 
an apparently rapid thermalization, like the existence of resonant Qq states [13] above 
Tc or values of transport coefficients close to the bounds predicted by the AdS/CFT 
correspondence U, 12, [3, [3, IH, [3], claimed to represent a good description of the 
sQGP for temperatures nearby Tc- See also Refs. 17, 18, 13, IS] for other approaches. 

At LHC on the other hand one expects to produce a fireball with a larger initial 
temperature (possibly Tq ^ 800 MeV) and a longer life-time. This could make weak- 
coupling calculations better justified. From the experimental side [2l|, the possibility of 
reconstructing also hadronic decays, like — > K~'k^ , in Pb-Pb collisions will remove 
the ambiguities encountered at RHIC. 

Two are the main issues addressed in this paper, both related to the transport prop- 
erties of charm and bottom quarks in a hot plasma. We first focus on the calculation of 
the transverse kt{p) and longitudinal hl{jp) momentum-diffusion coefficients. We then 
make use of the above results to follow the approach of heavy quarks to thermal equi- 
librium through a numerical Langevin simulation. We conffne ourselves to the case of a 
static medium, for temperatures spanning a range of experimental interest for LHC. 

In order to calculate the coefficients kt(2) and kl{p) we perform a generalization to 



the QCD case of an approach developed in [22, l23| to describe heavy (static) "quarks" in 



a hot plasma of electrons, positrons and photons. The latter represented in fact a system 
sharing important features with the QGP, but for which the possibility of employing in 
the calculations a Hard Thermal Loop (HTL) gaussian effective action allowed to study 
in an easy way important medium effects. This approach makes it possible to calculate 
also transport coefficients of QQ pairs. This is important, since measurements of pt 
spectra of quarkonia may allow to discriminate different scenarios. One expects in fact 
that the measured J ftp's in the final state should tend to follow the flow of the fireball, if 
they came from the recombination of charm quarks free of floating around in the plasma 
phase. On the contrary, if they remained bound and sufficiently close in space also in 
the deconfined environment, they would suff'er less collisions in the colored medium, as 
they would be seen mostly as a neutral object. 

As already mentioned, we then insert the above transport coefficients into a relativistic 
Langevin equation, allowing to follow the stochastic momentum evolution of the heavy 
quarks in the QGP. Their initial momentum po is taken of the order of the average 
transverse momentum p^p of c and b quarks produced in nucleus-nucleus collisions at 
RHIC and LHC. The latter can be, for instance, evaluated from a sample of events 
generated by PYTHIA [iS] with a proper tuning of the parameters. 

Dealing with the relativistic Langevin equation, with a momentum dependent diffu- 
sion term, is a non-trivial task. Actually, the latter belongs to the class of stochastic 
differential equations, which represents an open field of mathematical research, and re- 
quires a careful discretization procedure. In particular, the drag coefficient has to be 
appropriately tuned in order to lead, in the continuum limit, to a Fokker-Planck equa- 
tion admitting the relativistic Maxwell distribution as a stationary solution. We devote 
Appendix [Bl to discuss these issues in detail. 

Our paper is organized as follows. In Sec. Owe summarize well known results on the 
diffusion of non-relativistic heavy quarks in a hot plasma, displaying how the Langevin 
equation can be used as an effective theory to deal with the problem. However, in high- 
energy collisions charm and bottom spectra, though rapidly decreasing, have long power- 
law tails extending up to very large transverse momenta. Hence the need of extending 
the Langevin equation to the relativistic case. In Sec. Owe thus present the calculation of 
the transport coefficients for a quark with a generic momentum p. In Sec.|3]the approach 
is extended to deal with a QQ pair. Sec. [His devoted to the numerical results both for the 
transport coefficients and for the Langevin evolution of heavy-quark spectra. Finally, in 
Sec. m we discuss our results and draw our conclusions. Some technical details are given 
in the appendices. In Appendix |X] we recall the essential formulas for the HTL gluon 
propagators and spectral functions which are employed in the text. In Appendix |B] we 
give a self-contained discussion of the problems one has to face in solving the relativistic 
Langevin equation. In Appendix [C] we report the results of the kinetic calculation of the 
momentum-diffusion coefficient in the non-relativistic case, which is used as a benchmark 
to have some control on the uncertainties related to our effective HTL approach. 



2. Langevin approach for heavy quarks in the QGP: a brief summary 



The Langevin equation has been recently employed by several authors [25|, |26|, |27 



28l . I29L l3d | as an effective theory to describe the evolution in coordinate and momentum 
space of heavy quarks in a hot plasma. In this section we briefly summarize the main 
aspects of such an approach, focusing in particular on its range of validity and on the 
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information on the medium properties that can be extracted from the study of the heavy 
quark propagation. 

For the sake of simplicity our introductory discussion refers to the non-relativistic 
case, i.e. with the heavy quarks not too far from thermal equilibrium. This will be 
sufhcient for the above purposes. Actually, in heavy-ion collisions, c and b quarks are 
produced with a sizable transverse momentum, so that a relativistic study is in order. 
This however would introduce some technicalities; hence we postpone it to the following 
sections, where we shall see that, as long as T/E <C 1, the Langevin approach remains 
justified. 

In [i^jlllj it was shown that the interaction rate of a very massive (i.e. with M ^ T) 
quark at rest in a hot plasma is of order g^T and is given by (here we consider the QED 
case): 



r-g-r ^^ ,„, T... (1) 



Since the Debye mass mo gT is the only scale within the integral, the typical momen- 
tum exchanged in the interaction will be of its order0. In particular the maximum of the 
integrand occurs at q — tjid / so that most of the kicks received by the heavy quark 
will involve the exchange of soft electrostatic photons (since M is very large, one can 
neglect magnetic interactions and consider the static limit) with momentum, and virtu- 
ality, q ^ gT. Hence the duration of a collision Tcou ^ ^IgT will be well separated, at 
least in the weak coupling regime, from the average time interval between two scattering 
events At = 1/T ^ 1/g^T. Since Tcou <C At, one is then allowed to treat the diffusion 
of the heavy quarks as resulting from the sum of many uncorrelated momentum kicks. 

Due to the large mass Af, a huge number of collisions is required for an ensemble 
of heavy quarks to thermalize. At thermal equilibrium, from the equipartition theorem, 
one has Phcavy — v3MT, while for the light particles Piig^t ^ T , hence 

(Phcavy) "^(I'Ught)' (2) 

the same relation expressing the number of random collisions required to change the 
average squared-momentum by a factor of order one. In the case M ^ T the heavy- 
quark relaxation time will then be much larger than the one for the light particles of the 
medium, namely: 

M M 1 

Thoavy ~ ^ ^"s^t ^ g^T ln{l / g) ' 

where the last estimate refers to the weak-coupling regime (rught being of the order of 
the time interval between two hard collisions) and will be shown to hold in the following. 

In the Langevin approach [s^ the equation of motion for a non-relativistic heavy 
quark reads 

%--^DP'+eit), (4) 



Note that the total interaction rate is independent on the value of mo, which is only relevant to set 
the typical scale of exchanged momenta and to prevent infrared divergences by screening the interaction. 
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where the right hand side is given by the sum of a friction force (described by the drag 
coefficient rye) and a noise term, which is fixed by its temporal correlator: 



{Cm{t')) = K5^^5{t-t'). (5) 

The noise arises from the random uncorrelated momentum kicks received from the 
medium. Eq. ^ can be solved by discretizing the time derivative with a time-step 
Ai sufficiently large to include many collisions, but much shorter than the relaxation 
time of the heavy quark. The latter represents the typical time-scale over which one 
would like to follow the momentum evolution. The condition to be fulfilled is then: 

A /^N 
Tlight < At < Thcavy ~ flight- (6) 

The momentum diffusion coefficient k represents the averaged squared momentum ac- 
quired by the heavy quark per unit time, 

«=3(^), (7) 

and arises from the cumulated effect of many kicks suffered in the elementary time-step. 
It can be calculated from Eq. ([T]), by weighting the differential interaction rate with the 
squared momentum transfer, thus getting 

_ g'^mjjT q^dq 



67r Jo {q^ + mjj)'^ 
g'^m'jjT 



127r 



9max 9max 



(?2 



(8) 



Notice that, while the interaction rate is free of ultraviolet divergences — hard scat- 
terings not occurring very frequently — in the calculation of transport coefficients the 
latter play a major role, since, in spite of being quite rare, they can lead to a sizable 
momentum transfer. In the above gmax ^ T reflects an estimate of the maximum momen- 
tum exchange with a typical thermal particle. One often considers the result at Leading 
Logarithmic Accuracy (LLA): 

^LLA _ 5^"i|,T dq _ g'^mjjT 9,„ax 



^^^A ^ ^ -u- / n = ^ In:^, (9) 

67r J,nD 1 67r mo 

where the argument of the logarithm is of order 1/g, its precise numerical value repre- 
senting the theoretical uncertainty of the present approach. 

By requiring the momentum distribution to reach thermal equilibrium one gets the 
Einstein relation between the drag and momentum diffusion coefficients: 



^ 2MT- 



Concerning the evolution in space, from 



{x^{t)) - 6Dt, with x'{t) = / dt'i^^, (11) 

t^oo _/„ M 
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one derives the (spatial) heavy- quark diffusion coefficient 



D = -!—^—, (12) 
MrjD K 



whose expression at LLA follows from Eq. ^ 

VlixT 



jjLLA ^ . (13) 



niD 

For the heavy-quark relaxation time one has then: 



Thoavy = — - ^D. (14) 

By comparing Eqs. Q and (HU one gets the interesting result [j^t 

s«i 54rin(l/g) 

showing that the heavy quark diffusion coefhcient gives a reasonable estimate of the 
relaxation time of the medium. We also notice that, in the weak coupling regime, from 
Eq. in]), one has rught ~ l//Tln(l/g). 

For a recent study of the possibility of estimating heavy-quark transport coefficients 
from euclidean lattice simulations see Ref. [sst . 

3. Propagation of a heavy quark in a hot plasma: transport properties 

We start our investigation by considering the propagation of a heavy quark in a hot 
ultra-relativistic plasma. When the quark is very massive and/or is endowed with a very 
large momentum p ^ T, it is reasonable to describe it within the eikonal approximation, 
in which it moves along a straight-line trajectory, acquiring a phase due to the interaction 
with the background gauge field. 

One could argue that, due to the multiple kicks received from the particles of the 
medium, during its propagation the heavy quark would loose energy and acquire more 
and more transverse momentum, so that at some point the assumption of straight-line 
propagation should cease to be valid. However, one can take advantage of the huge 
separation, occurring for Ep ^ T, between the relaxation time of the medium, rught, and 
the much larger time required by the heavy quark to approach equilibrium. In order to 
evaluate transport coefficients (related for instance to heavy-quark energy loss, diffusion 
and momentum broadening), it is sufficient to follow the propagation of the quark for 
a time large compared to riight but still much shorter than the relaxation time of the 
high-momentum heavy quark. For this purpose our eikonal approach should be justified 
and indeed turns out to provide results in agreement (at least at Leading Logarithmic 
Accuracy) with the ones obtained in other approaches as, for example, by solving the 
Maxwell equations in a dielectric medium [sj, [s^ [s^. If one is really interested in 
following the relaxation of the heavy quarks toward thermal equilibrium, then one can 
use the above findings for the transport coefficients and insert them into the Langevin 
or Fokker-Planck equations. 
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Following the approach developed in Refs. 22, l23|, the configurations of the gauge- 
field entering into the eikonal phase are weighted by the HTL effective action. We start 
from the QED case (i.e. a plasma of photons, electrons and positrons), in which the 
latter is gaussian; this allows to perform the functional integral exactly, leading to an 
exact exponentiation of the gauge-field propagator. Later we will show how the results 
for the transport properties can be generalized to the QCD case, yet recovering well 
known findings. 

In analogy with the approach adopted in Refs. 37 1, we consider the retarded 
propagator of a heavy quark — which is treated as a test particle — created at {0,r[) 
and annihilated at (t,ri); it is defined as 

G^(i,ri|0,r;) = t0{t){ij{t,r,)^\O,r[)), (16) 

where the expectation value refers to a thermal average over the states of a hot medium 
of light particles. In the eikonal approximation it is given by: 

G^(t,ri|0,r'i) ^i0{t)6{ri-r[-vt)G{t). (17) 

In the above, the eikonal phase (which is complex, the imaginary part accounting for the 
effects of collisions) is expressed in terms of the real-time gauge-field propagator given in 
Appendix El 

' ''d^x [ d^yJ'^{x)D^,{x-y)r{y) 



G{t) = exp 



(18) 



2 

where the current describing the propagation of the heavy quark is given by: 

J^(a;) = 90{x°)9{t - x°)6{x ^ r[ ^ vx°){l,v). (19) 
After expressing the gauge propagator in Fourier space one gets: 

— r i 2 f duj f 2 [1 - cos(w - q ■ v)t] 
G{t)=e.p{-g J^Jj^ (u^-q-vr 

X [Dl{uj, q) +v^l-{v- qf) Dt{uj, q)] }, (20) 

whose large-time behavior can be obtained from the limit 

1 — cos(ijj — q ■ v)t N /„ X 

lim ^ ^i--^ = 7rM(cj - q • t> . (21) 

t-+oo [lj — q ■ vy 

In particular the probability of finding the heavy quark with momentum p — ^Mv will 
tend to decrease due to the collisions with the plasma particles, whose effects are encoded 
into the imaginary part of the HTL gauge-field propagator (see Appendix |^ : 



ImDi/j^(w,(7) = pL/T{uJ,q) (^{uj) -f ^ 



(22) 



This allows to formally define an interaction rate for a quark propagating with velocity 
V (taken in the following along the z— axis), which turns out to be given by (the spectral 



functions being odd) 



duj 



duj 



dq 



(27r)3 
dq 



6{lu 



v) [pLito, q) +v'{l-{v- qf) priiu, q)] N{uj) 



q ■ v)p{uj,q)N{uj) 



(23) 



Note that, strictly speaking, the above expression develops an infrared divergence aris- 
ing from the exchange of long-wavelength magneto-static gluonH. However the above 
processes cannot affect the transport properties related to the heavy-quark propaga- 
tion (energy loss, momentum broadening, ...), since they are related to negligible en- 
ergy/momentum exchanges. Hence we can take advantage of Eq. (j23p for the computa- 
tion of the transport coefficients we are interested in. 
The energy loss per unit length is then given by: 



dE 
dx 



V 



duj 



dq , 
J^^5{uj - q ■ v)p{uj, q) ujN{uj). 



(24) 



S{uj — q ■ v) 



(25) 



The Dirac delta 

— si cos 9 — — 

qv \ qv^ 

can be exploited to perform the angular integration. Furthermore, since the spectral 
bmction is odd, one can replace the Bose distribution with its even part 



N{uj) 



N{lu)+N{-lu) 



1 

2' 



(26) 



thus obtaining 

dE 
dx 



dqq 
qdq 



duj i 



PL{uJ,q) + V 



PT{oJ,q) 



duj ujp{uj, g). 



(27) 



where, in the high-energy limit, 



--^ VET, with E the energy of the heavy quark, 
is the maximum momentum transfer in a collision with a plasma particle, whose typical 
momentum is of order T. At LLA the energy loss can be computed analytically, by 
setting a lower bound of order m/^ in the momentum integration and using for the 
spectral functions the expressions given in Eq. (j76p . One gets: 



dE 

dx 



LLA 



9 9 

g mf) 

9 9 



mo 

9ma 



dq 
dq 

q 



dio LO 



dxx 



UJ^/q^ 



2{1- uj-^/q^) 

„,2 J2 -\ 



2(1 -x2) 



(28) 



■^Actually, lattice QCD results provide evidence for the existence of a non vanishing magnetic 
mass [SSI . of non-perturbative origin, which would eliminate the above divergence. This in principle 
could affect our results for the transport coefficients, but such an issue lies beyond the scopes of the 
present paper. 
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After performing the integrals one gets: 



dE 
dx 



g^m^D gmax 



LLA 



■niD 



1 



1 - i;2 1 

In — 

2v 1 



(29) 



to be compared with the analogous result displayed in Eg. (B 31) of Ref. [26| for the QCD 
case and quoting previous findings obtained in Refs. [39j, |40|. 

We now address the calculation of the momentum diffusion coefficients kt and kl 
which will enter the noise term in the Langevin equation [see Eqs. (jlOip and (jl02p ]. 

Let us start with the transverse momentum diffusion coefficient, representing the 
mean squared transverse momentum acquired per unit time by the heavy quark crossing 
the medium, namely: 

1 Ap|, 
2 



At 



(30) 



Also the latter can be computed starting from Eq, 
■-2 ^ r dq 



9_ 
2 



duj 



S{uj-q- v)p{uj, q)q\l - cos^ 6')7V(w). 



(31) 



In the above only the odd part of the Bose distribution contributes to the integral, so 
that one can make the replacement 



iV(cj) 



7V(w) 



— ^ = -coth--, 



(32) 



thus obtaining: 



dqq \ dujp{uj,q)q^ 1 



coth 



(33) 



In analogy with the energy loss calculation, the LLA result for the transverse momentum 
diffusion is obtained by employing the corresponding approximate expression for the 
HTL spectral functions given in Appendix [XI setting an infrared cutoff of order m^i in 
the momentum integration and assuming that the relevant contribution to the above 
expression is given by processes with small energy transfer, so that one can approximate 



coth 



(3uj 



2T 



(34) 



One obtains: 



^LLA 



7 7 

g^Tra] 



dq 



dq 

q 



dlO LJ 



dx 



uj^/q 



2 1 



2(l-tj2/9^) 



1 



2T 



2(1 -a;2) 



After performing the integrals one gets 



,LLA_ 



g^Tm 



D 



Ait 



In- 



1 

2^ 



(35) 



(36) 



again in agreement with the leading log term in Eq. (B32) of Ref. [26|. Note that 
the transverse momentum diffusion is closely related to the transport coefficient q, often 
entering into the radiative energy loss calculation as a parameter encoding the properties 
of the medium, and representing the transverse squared momentum acquired by the 
propagating parton per mean free path 41|. One can in fact measure the distance X 
covered in the medium in units of Amfp, setting X = nAmfp, so that 



AX = X 



Hence 



q ■ 



_ {Ap: 



,Ap 



infp ■ 
V 



Amfp AX v' At 
leading, in the case of an ultra-relativistic parton, to 



-KT, 



g'^Tmj 
2tt 



In 



mo 



in agreement with the findings of Ref. [42| • 

The calculation of the longitudinal momentum diffusion 



^ At 



follows in perfect analogy. From 



one gets 



KL= 9 



dw 



dq 



S{ll> - q ■ v)p{uj, q) q^ cos^ N{u)) 



qdq 



vq 



duj p[uj, q)^ coth . 

v-^ 2 



Again the result at LLA can be expressed analytically: 



.LLA . 



dq 
7a 



g^Tm 



mo 



dq 

q 



—ri dhJOJ 



2/2 



dx X 



2(1- 0-2/92) 
0,2 ^2 n 



2(1 -a;2) 



which leads to: 



,LLA 



frmi^^q 
47ru2 



ruD 



1 



1 



■In- 



1 



(37) 
(38) 

(39) 



(40) 



(41) 



(42) 



(43) 



2v -1-..' (44) 

to be compared with Eq. (B33) of Ref. [2^1 ■ Note that at LLA the relation between 
friction and momentum-diffusion coefficients, required by the Langevin approach (see 
Sec. Hand App.[B|), holds: 



,LLA 



2T 


dp 


2T 


dE 




V 


'dt 


V 


dx 


LLA 



2TErj^^, 



(45) 



as it can be seen by comparing Eqs. P5|) and (UH). We notice however that, as it will be 
discussed in the following, when the strength of the noise depends on the momentum of 
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the brownian particle (as it turns out to be the case) the viscous term in the Langevin 
description needs a correction (subleading in T/E) in order to recover the correct con- 
tinuum Umit to the Fokker-Planck equation. 

Let us now check a posteriori that in a time interval of order riight the momentum 
loss Ap and the acquired y/Jj^ are negligible with respect to the initial Ep ^ T, so that 
the calculation of transport coefhcients within the eikonal approach results meaningful. 

From dp/dt = dE/dx, during the time rught, one gets: 



|Ap| = 



dE 




dE 


dx 


■''light ^ 


dx 



D 



g'^mjj gmax 



UttT 



8tt 



T-^p. 



Concerning the transverse momentum, during the time interval rught, one has 

(Apl) = 2KTTiight ^ ktD, 



(46) 



(47) 



so that 



47r 

V 

T<^p. 



1/2 



In^ 



U-kT 



^max 

mo 



(48) 



4. Transport properties for a QQ pair 

Within the same framework it is possible to study the transport properties of a 
quarkonium, which we simply model as a QQ pair, of a fixed size, propagating with 
velocity d in a hot plasma. Similar studies, limited to the ener gy- loss problem and to 



the definition of a "dipole potential", can be found in Refs. [3J, |36|. Here the problem 
is addressed by solving the Maxwell equations in the linear response approximation for 
a QQ pair (playing the role of an external current) propagating in the QGP, modeled as 
a dielectric medium. 

Clearly, the results will depend not only on the size of the "dipole" , but also on its 
orientation with respect to the direction of propagation. Following the choice adopted 
in Rcf. [3(|, we take, without loss of generality, v along the 2;-axis, r= {rj_,0,z) in the 
2:a;-plane and the exchanged momentum as <j' = (7(sin^ cos sin0 sin0, cos^). 

We consider then the creation of a quark at (0, r[) and of an anti-quark at (0, ra), with 
r = r'l — r'2, and we follow the propagation of this pair within the eikonal approximation, 
as already done for the single particle case in Sec. [H Hence, one can take advantage of 
Eq. (jlSp. where now the current describes the propagation of a dipole and is given by, 

Jf (x) = ge{x°)9{t - x°) [S{x -r[- vx") - 6{x - - -ux")] (1, v), (49) 

so that for the QQ propagator one gets: 

duj f dq 2[(1 — cos(ijj — q ■ v)t] 



G{t) = exp|i,g 



27r J (27r)3 {uj - q ■ v)'^ 
X [1 - cos(q • r)] [Dl{u;, q) + v' {l ~ (v ■ q)^) Dt{uo, q)] }. (50) 
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One can again consider the large-time decay of the above correlator, which allows to 
identify the interaction rate of the dipole: 



r = 2.g" 



dq 

(2^ 



[1 — cos(q • r)] 5{ijj — q ■ v)p{iLj, q)N{uj). 



(51) 



The Dirac delta can be used to perform the integral over the polar angle, while the 
azimuthal integration gives rise to a Bessel function, by exploiting the integral represen- 
tation ^ 

d(f> j.^ coscj) 
27r 



(52) 



with q± — qsinO. One can then study the transport properties for the QQ pair. 
We first examine the energy loss which is given by [see Eq. ([27]) ] : 



dE 



9 



qdq 



du 



1 - Jo qrj 



'1 



q2y2 



dx 2t:^v^ 

For the transverse and longitudinal momentum diffusion coefficients one gets 



p{uj,q)uj. (53) 



9 



qdq 



x'p{uj, q) coth 



duj 



(3uj 



I - Jo ^qrj_^l 

1 



q2y2 



q2y2 



(54) 



and 
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qdq 



xj){uj, q) coth 



duj 



l3uj\ LO 



^ - Jo \ qr l 



q2y2 



(55) 



respectively. 

In the calculation one also needs to account for the dependence on the dipole orien- 
tation. For this purpose we start by considering the rest frame of the propagating pair 
(whose coordinates we label with a tilde): for a 5— wave state (till very large velocities 
the in-medium potential is still approximately spherically symmetric [43]) one has 



Moving then to the rest frame of the thermal bath (the laboratory frame), if the pair 
propagates along the z— axis, its transverse size (in the xy— plane) is not affected and 

one has: 

r^ = ^{x^) + (y2) = + (y2) = y^i?. 

On the other hand, the longitudinal size of the dipole turns out to be Lorentz-contracted, 
according to the relation 

z = S/j^ R/{V3j). (56) 
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The above expressions of rj_ and z will be employed in the calculation of the transport 
coefficients. 

The small dipole limit can be easily estimated by considering the expansion 
I — Jo{qri_fimO)coB{qzcoBO) ^ — 1 

= ^(^sm^e + ^j; (57) 

hence, for small sizes, the transport coefficients grow linearly with the transverse area of 
the dipole. One gets then, in the limits of very small and very large velocities, 

dE\ g^R^ r . N 

^ ' q^dql dcjp(tJ,q) 1-— w (58) 



for the energy loss, 

^''^^^«' H«;-L "''^"l ^'^^(^'«)^°th( — )g2(i- 

(«t)„^1 fj^\'dqj' cL:-p{u,q)coi\,{^)q'{l-^\ (59) 

and 



^2^2 /-S 

~ -TTTT' / 9 "9 / at^p(t^, (?) coth ' ' ■ 



for the transverse and longitudinal momentum diffusion coefficients, respectively. Notice 
that for V —^ 1 terms containing z are neglected, since they are suppressed by the Lorenz 7 
factor: in the ultra-relativistic limit only the transverse area is relevant for the transport 
properties. In practice, however, the above expansions are not so useful, since the results 
depend quadratically on the momentum cutoff g,nax, which is not well determined and 
gives a huge systematic theoretical uncertainty. On the contrary working with the exact 
formulas provides transport coefficients with only a mild logarithmic dependence on the 
ultraviolet cutoff. 



5. Numerical results 

Here we present our numerical results, starting from the transport coefficients kt{p) 
and kl(jp) introduced in Sees. [3]and[ll they are then employed in the Langevin evolution 
of the momenta of a large sample of heavy quarks (charm and bottom). The rigorous 
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results obtained in Sec. [3] are generalized to the QCD case by adding a color charge to 
the current given in Eq. (I19p . setting 



JA'«(a;) ^ q''g6{x^)e{t - x^)5{x - ri - vx^){\,v), (61) 

with a = 1, . . . ; iV^ — 1; the corresponding transport coefficients can be simply obtained 
by multiplying the ones derived in Sec. [3] by the Casimir factor Cp = In the 

present exploratory study the QGP is modeled as a uniform static medium, considered 
at different temperatures spanning a range of interest for heavy ion collision experiments. 

As already mentioned, the results obtained within the present framework display a 
logarithmic dependence on the ultraviolet cutoff gmax, related to the maximum momen- 
tum exchanged in a collision with a typical thermal particle. Actually, this represents an 
intrinsic uncertainty of our scheme, which — in dealing with large momentum transfer 
processes — should be supplemented by a microscopic kinetic calculation involving the 
evaluation of the Born matrix elements for the scattering of the heavy quarks with the 
particles of the medium However, for the sake of ease and consistency we 



prefer to pursue the approach developed in Sec. [S] indeed, it is able to catch all the 
essential qualitative features of the transport coefficients (e.g. the growth of k^/l with 
the momentum p, its suppression in the quarkonium case, etc.) and is simple enough 
(also in view of a more realistic study of an expanding medium) to be implemented into 
a Langevin simulation. Bearing in mind that the ambiguity related to the choice of gmax 
is unavoidable, we refer the reader to Appendix [Cl for details on the procedure employed 
by us to fix a realistic value for it. 

Concerning the value of as, which one expects to be quite large at RHIC and to 
approach a weak-coupling regime at the (larger) LHC temperatureSjit has been taken 
from the 2-loop QCD /3- function with the parameters given in Ref. f44| . One has 

g-\^i) = 2 6o In (m/Aqcd) + ^ In [2 ln(A*/AQCD)] , (62) 

with 

- 16^ - ' - (16^ - f ^0 (63) 

and Aqcd = 261 MeV. The QCD coupling was then evaluated at a scale proportional to 
the temperature. As two representative cases we present the results corresponding to the 
typical choices /i = ttT and /i = 27rT, the differences reflecting the lack of predictivity 
of the calculation. A more careful study of the running coupling effects (actually in the 
evaluation of coUisional energy- loss) was given in Refs. [4^13], but this goes beyond the 
scope of the present analysis. 

We also compare our results, both for the transport coefficients and for the Langevin 
dynamics of the heavy quarks, with the findings provided by the AdS/CFT correspon- 



dence, which, for the A/" = 4 SYM theory, predicts: [l3|, [ij, ll5| 



Kt = 7 ' VATTJgYM, 

KL = t'/'VA^T^Iym, (64) 



where 7 = l/\/l — v"^ and A = 5gYM^c- Actually, translating the above results to the 
hot-QCD case is not free of ambiguities and different strategies were proposed. Here we 

14 



HTL (n=2jt) 
LLA (h=2ji) 
HTL (ii=7i) 
AdS/CFT 




T (GeV) 



Figure 1: The HTL momentum-diffusion coefficients Kj^/l divided by T'' for a bottom quark of mass 
M = 4.2 GeV and momentum p = 4 GeV as a function of the temperature. The HTL results, for fi = 2it 
and n = TT, are shown (solid and dashed lines, respectively). The dot-dashed lines correspond to the 
LLA results. Thick lines refer to kj-, thin lines to k^. The AdS/CFT values are also reported. 



follow the recipe adopted in Ref. , which is based on matching the energy density and 
the QQ force in the two theories, namely: 

S^/^TsYM = Tqcd and A = 5.5, (65) 

the latter being chosen within the range 3.5 < A < 8 allowed by the above mentioned 
procedure. It was pointed out in Ref. that the results in Eq. only hold for 
heavy-quark momenta corresponding to 7 factors below a critical limit: 

\ V AJsYM/ 

To our knowledge no calculation is available for larger momenta. Hence, in our plots, 
as 7 exceeds this critical value, we freeze kt and kl to their estimates below jc- Had 
we trivially identified the coupling and the temperature in the two theories, requiring 
A = Gtt to reproduce a QCD plasma with Nc = 3 and as = 5^/47r « 0.5, the range of 
validity of the AdS/CFT calculation would have been even more limited. 

In the range of temperatures covered by our analysis the transport coefficients k-^ 
change dramatically, increasing approximately as (deviations being due to running- 
coupling effects) , as it is displayed in Fig. [TJ This leads to observable consequences for 
the heavy-quark dynamics, that we will discuss. Fig. [1] refers to the case of a 6-quark 
of not too large momentum, so that the condition in Eq. ([M]) is satisfied. In the same 
figure we also show for comparison the LLA result for /i = 27r. 

In Fig. [2] we give, for a charm quark and for different temperatures, the coefficients 
kt{p) and i^l{p) which are derived from our HTL calculation. We also display, for 
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Figure 2: The HTL momentum-diffusion coefficients k.t/l{p) for ^ charm quark of mass M = 1.2 GeV 
for a range of temperatures above Tc of experimental interest. We also plot, in the regime where they 
are available, the AdS/CFT results. Again, thick lines refer to kt and thin ones to k^. 




2 4 6 
p (GeV) 



4 6 8 
p (GeV) 



Figure 3: As in Fig. [S] but for a bottom quark of mass M = 4.2 GeV. 
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Figure 4: Left panel: k,t,l{p)/ t^T,L{0) as a function of momentum for charm quarks (rric = 1.5 GeV) 
at T = 200 MeV in our model (HTL), in the the QCD kinetic calculation of Ref. '2(| (LO QCD) and in 
the resonance model of Ref. ^10|. Right panel: the drag coefficient rjoip) as a function of momentum for 
charm quarks (mc = 1-5 GeV) at T = 300 MeV in our model (HTL), in the calculations of Refs. fT3.[20ll 
(pQCD) and ^ (LO QCD). 



comparison, the corresponding AdS/CFT result. As expected, the latter model provides 
stronger coefficients and also a much steeper momentum dependence of k,l{p)- Notice 
that, for large enough temperatures, 7c < 1, so that the condition 7 < 7c cannot be 
satisfied. In Fig. [3] the same quantities for the case of a bottom quark are shown. 

In Fig. m we compare our results to a few of the calculations for the transport coeffi- 
cients available in the literature. In the left panel we display the ratios kt.t.Ip) / kt. l(O) 
at r = 200 MeV in our HTL model together with the QCD kinetic calculation of Ref. [26| 
(LO QCD) and the resonance model of Ref. [l3|. The momentum dependence of the lon- 
gitudinal coefficient has on the whole a similar trend in the HTL and LO QCD models, 
whereas in the latter the transverse coefficient displays a less pronounced momentum 
dependence. On the other hand, the resonance model provides coefficients that are es- 
sentially flat in the available range of momenta. 

In the right panel of Fig. |4] we compare the drag coefficient 77^1 as obtained in the 
HTL model at T = 300 MeV to the LO QCD model of Ref. '2(il- In that model the value 
of i]d{0) is parameterized through the diffusion constant in space D, which is treated as 
a free parameter, and in the figure we display two choices for D in the range considered 
in Ref. . As one can see, the momentum dependence is similar to the one of the HTL 
model, the absolute normalization covering a wider range of values. 

In the same panel we also report the results obtained in Refs. where an 

emended perturbative QCD (pQCD ) calculation is employed by introducing a smaller 
infrared regulator and a suitably chosen running coupling constant. In the figure the 
labels (C) and (E) refer to two different choices of parameters (see Table I of Ref. [S])- 
The coefficient r]u is smoothly decreasing with p in this model, whereas in the HTL 
model (and also in the LO QCD one) it displays a maximum. At momenta p > 3 GeV 
the magnitude of rju is similar in the HTL and pQCD calculations, the latter, on the other 
hand, being definetily larger at small momenta. Also the resonance model of Ref. 
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R (fm) p (GeV) 

Figure 5: The charmonium transport coefficients, normalized to the result for two independent c-quarks. 
Left panel: as a function of the QQ separation, for a pair propagating with momentum pcc = 2pc = 2 
GeV. Right panel: as a function of the quark momentum, for a pair of size R = 0.4 fm. Thick lines refer 
to Ky, thin lines to ft^. Various temperatures are considered. 



displays a smoothly decreasing behavior with momentum. 

We then consider, in Figs. [5] and [H the case of a QQ pair crossing the QGP. The 
quarkonium is modeled as a pair of heavy quarks, which propagates in the medium with 
a fixed separation and with a momentum given by the sum of the two individual quarks. 
Any binding effect is neglected, hence one has simply 

^ {El^pl (67) 

with AU = 2Mf w 2 • 1.5 GeV in the case of J/V" and = 2M°^ w 2 • 4.7 GeV in the 
case of a T meson. 

In spite of its simplicity such a model is able to display how — due to the destructive 
interference effect encoded in the [1 — cos(q • r)] factor of Eq. ([5T|) — quarkonium states 
of small size propagating in the QGP would suffer much less rescatterings with respect 
to two uncorrelated heavy quarks, as it appears from the left panels of Figs. [5] and [H In 
particular, since the relevant quantity is the ratio between the average separation R of 
the heavy quarks and the Debye radius Rd = for a given dipole size the suppression 
is more dramatic at small temperatures, where the Debye radius is larger, implying a 
smaller value of R/Rd. 

In the right panels wc still consider the quarkonium transport coefficients, but as a 
function of the momentum p of the single "constituent" heavy quark. We take R = 0.4 
fm in Fig. [H]and R = 0.2 fm in Fig. [HI corresponding to the typical mean square radius 
of the charmonium and bottomonium ground states, respectively. 

We now make use of the above findings for kt/l{p) to follow the relaxation to thermal 
equilibrium of a large sample of charm or bottom quarks. We model such a process with 
the following Langevin equation: 

'^^-VD{p)f + C{t). (68) 




R (fm) p (GeV) 

Figure 6: The bottomonium transport coefBcients, normalized to the result for two independent 6-quarks. 
Left panel; as a function of the QQ separation, for a pair propagating with momentum p^;; = 2pi, = 2 
GeV. Right panel: as a function of the quark momentum, for a pair of size R = 0.2 fm. Thick lines refer 
to kt, thin lines to /t^. Various temperatures are considered. 



As in the non-relativistic case introduced in Sec. [21 the resulting stochastic evolution of 
the heavy-quark momenta is determined by the noise correlation function 

{C{t)e{t')) = Sit - 1') [KL{p)fp' + KTipw - fpn] , (69) 

which is expressed in terms of the transverse and longitudinal transport coefficients pre- 
viously evaluated. In Eq. ((55)) the drag coefficients r^nip) is fixed in order to ensure the 
approach to equilibrium. An extended discussion of the relativistic Langevin equation, 
of its discretization and of the algorithm employed in its numerical implementation is 
reported in Appendix [Bl We mention here that other choices can be found in the lit- 
erature. Some authors, for instance, performed a first-principle calculation for rjD [29| . 
fixing the remaining coefficients by requiring the proper equilibrium limit. Our choice 
is motivated both by the quest for consistency (since the same microscopic calculation 
provides both kt and kl) and by the importance of distinguishing between the trans- 
verse and longitudinal momentum broadening of a fast particle (at variance with what 
done in [1^). Finally, we observe that, in any case, in dealing with a stochastic differ- 
ential equation, the - momentum-dependent - friction coefficient receives a correction 
depending on the adopted discretization scheme, so that it appears more natural to fix 
the latter "by hand". We point out, however, that the HTL calculation, at least at LLA, 
provides consistent results for rio and kl- 

We start with the case of a sample of 2 • 10® charm quarks, with an initial mo- 
mentum po representative of their typical initial pr in nucleus- nucleus collisions. The 
parametrization of their spectrum used in Ref. [26|, referring to RHIC conditions, gives, 
for instance, prp = 1.23 GeV. In the case of Pb-Pb collisions at LHC (at ^/smn — 5.5 



TeV), one gets, instead, pj- = 2.18 GeV [2l| for charm quarks generated by PYTHIA 



with parameters tuned to reproduce the results by Mangano et al. [47| ■ The momentum 
distributions resulting from the Langevin evolution, for different values of T and po, are 
presented in Figs. [7] and [51 for c-quarks at increasing values of time. For any value of 
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T, for large enough times, the charm momenta turn out to be described by a relativistic 
Maxwell- Jiittner distribution 



JmAp) = n 7 r-r, with / (fpfuTip) = 1, (70) 

where K2 is a modified Bessel function. We notice, however, that the approach to 
equilibrium is much faster for larger values of the temperature, due to the huge increase 
of the coefficients kt/l{p) with T displayed in Fig. [T] 

For comparison we also give, in the right panels, the corresponding results arising 
from the AdS/CFT estimate of the transport coefficients entering into the Langevin 
equation. Such a strongly-coupled scenario gives rise to an extremely fast relaxation 
toward equilibrium. Notice that the case T — 800 MeV is out of the range of validity of 
the present available methods to compute kt/l{p) within the SYM framework. 

We then consider, for the same range of temperatures, the Langevin dynamics of an 
ensemble of 6-quarks with mass M = 4.2 GeV, initialized with a momentum po = 6 GeV, 
chosen of the order of the average px of the spectrum generated by PYTHIA, with the 



parameters given in Ref. [21|. Notably this corresponds to a much harder momentum 
than the equilibrium value. This fact, together with the larger mass of the bottom quark, 
makes the approach to equilibrium much slower with respect to the case of charm. Our 
results for the momentum distribution function are displayed in Fig. [51 As it can be seen 
the HTL calculation provides relaxation times too large to be of interest for the realistic 
case of an expanding (and cooling) fireball. 

On the contrary the AdS/CFT scenario supports a quite rapid thermalization also 
for 6-quarks, the latter being almost immediate at the largest temperatures examined. 

Finally, in order to get a feeling of the rapidity of the thermalization process in the 
different cases considered here (which can be not so easy to grasp just by looking at the 
distributions), we give in Figs . fTOl and 1 11 1 the evolution of the mean squared momentum 
(p^) to its equilibrium value predicted by the Maxwell- Jiittner distribution of Eq. (fZO)) . 
The difference between the two models (HTL and SYM) appears particularly evident in 
the case of 6-quarks, which could be an interesting probe to study at LHC in order to 
discriminate between different scenarios. 



6. Discussion and conclusions 

The two main issues addressed in this paper are the evaluation of transport coefficients 
for heavy quarks in the QGP and the study of the resulting Langevin evolution of their 
momenta. We showed that important qualitative and quantitative (though affected by 
some systematic uncertainty) information on the energy-loss and momentum-broadening 
of heavy quarks can be extracted from an approach (based on the HTL approximation) 
developed for a different purpose, namely for the definition of an effective in-medium 



potential for a pair of static quarks [2^, |23| . 

Concerning the transport coefficients, both kt and kl were found to display a sizable 
growth with the momentum, increasing faster than kt- 

The approach developed in this work has been also useful for the study of transport 
properties of quarkonia. QQ states of small size, behaving as an almost neutral object, 
turned out to suffer less rescatterings. This could be of phenomenological interest in 
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Figure 7: The momentum distribution at diflferent times t of charm quarks with M = 1.2 GeV in a QGP 
for a range of temperatures. The left and central panels display the HTL results {/M = 2'kT and [j, = ttT, 
respectively), while the right panels refer to the AdS/CFT calculation. For the initial momentum we 
take po = 1 GeV. Square-dots represent the Maxwell-Jiittner equilibrium distributions. 
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Figure 8: As in Fig. [T] but for an initial momentum po = 2 GeV. 
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Figure 9: As in Fig. [7] but for a sample of bottom quarks of mass M = 4.2 GeV and initial momentum 
po = 6 GeV. 
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Figure 10: The approach of (p^) to the equihbrium value predicted by the Maxwell- Jiittner distribution 
in the case of charm quarks at various temperatures. 
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order to discriminate between different production/suppression scenarios. Analysis of 
momentum spectra would, in fact, help to understand whether the observed quarkonia 
in the final state arise from the recombination of uncorrelated heavy quarks or from their 
primordial production, after crossing unaffected the deconfincd environment. 

We then addressed the Langevin evolution of a large sample of heavy quarks (charm 
and bottom) — starting from an initial momentum chosen to be representative of the typ- 
ical ones in nucleus-nucleus colliwions — and following their approach to thermalization. 
We also tried to give some basic discussion of the essential aspects of the relativistic 
Langevin equation as an effective model to describe the relaxation of an ensemble of 
heavy particles toward equilibrium. 

Our numerical findings allow to draw some general considerations. 

We first notice that at the largest temperatures explored (T = 600 and T = 800 
MeV) the thermalization of c-quarks predicted by our (weak-coupling) HTL calculation 
is quite fast. Such an occurrence could be of phenomenological interest for LHC, where 
one expects quite large temperatures to be achieved in the first instants after the collision. 
On the contrary the relaxation of 6-quark spectra toward equilibrium arising from the 
present HTL calculation is always too slow to allow them to follow the flow of the medium 
in the realistic case of an expanding fireball. 

For what concerns the strongly-coupled AdS/CFT scenario, c-quarks are found to 
thermalize almost immediately, in spite of our conservative choice of freezing the values 
of K-T/L for 7 < 7c. Interestingly, also 6-quarks approach equilibrium quite fast (at the 
largest temperatures), making possible for them to inherit at least part of the flow of 
the medium produced in heavy-ion collisions. Elliptic-flow is in fact expected to develop 
in the very initial stage of the fireball evolution, when the medium is extremely hot, 
entailing large values of Kt/l (which, we remind, grow as T^). 

At the large temperatures realized at LHC, &-quarks appear then as a promising 
probe to discriminate between the weakly and strongly-coupled scenarios. A sizable flow 
of 6-quarks could not be compatible with any weak-coupling perturbative calculation. 
Measurement of the elliptic-flow and of the quenching of pr spectra of single-electrons 
(or positrons) arising from the semi-leptonic decays of B-mesons (which will be possible 
at LHC) could then shed light on the transport properties of the QGP at the future 
accessible regime of temperatures. 

In future work we plan to employ the approach developed here, whose numerical 
implementation turns out to be quite inexpensive, in order to address the realistic case 
of a fireball displaying longitudinal and transverse expansion. Supplementing our study 
with a more microscopic kinetic calculation of the hard-scattering contribution to the 
transport coefficients is also object of our current investigation. 
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A. The HTL propagators and spectral functions 

We give here the explicit expressions for the HTL gluon (photon in QED) propaga- 
tors in the Coulomb gauge, together with their spectral functions, which are employed 
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throughout the text. Further details can be found, e.g., m Ref. |48!|. 

The longitudinal and transverse analytical (off the real-energy axis) propagators read, 
respectively, 



At(9°,'z) 



-1 



-1 



(71) 



where 



nUx) = ml[l^Q{x)], 

? 

TIt{x) = [a;^ + (1 - x'^)Q{x)] and 

, . _ X x+1 



(72) 



being x = q'^ /q. 

The corresponding retarded propagators are obtained by setting q'^ = to + iri in the above 
expressions, namely: 

D^^^{LU,q) = AL/TiL^ + tV,q)- (73) 
The HTL gluon (photon) spectral function follows immediately from the definition: 



PL/ri^.q) = 2ImL»f/-r(tj,q). 
One gets (for space-like momenta): 



(74) 



PL{uj,q) 



PT{^,q) = 



2q 



UJ + q 



LU — q 



'2q 



2 Oj{(f-<J^ 
2g3 



(75) 



2 2 

— q 



l-^^ln 



2u}q 



UJ — q 



LjjiJ^ - q^) 
' 2g3 



Their expressions at LLA are obtained by neglecting the self-energy corrections in the 
denominators and are given by: 



PLi(^,q) 



2 UJ 

ZTTrrLn — 
^2(7 



LLA 



nmj 



PT{uJ,q) 



2q 



LLA q^{q'^-LJ^)' 



(76) 



This implies that the infrared divergences are no longer screened by the Debye mass, and 
must be regulated by hand by setting a lower cutoff of order mjj in the integration over 
the three-momentum of the exchanged gluons (photons). 
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When dealing with the in-medium heavy quark propagation, within the eikonal ap- 
proximation, one needs to consider the exponentiation of the real-time gauge-field prop- 
agator; the latter is given by 

DL(i^,q) = ^ + / 5 —— +ipL[uJ,q)N{uj) 

q J-oo 27r - (w + iVi) 

DT{uJ,q) = / 5 -—— + tpT{uJ,q)N{uj), (77) 

where N{uj) is the customary Bose distribution. The above expressions are related to 
the respective retarded propagators as follows: 

Dl/t{^, q) = ^f/T(^> q) + *Pl/t(w, q)Niij) (78) 

= ReD'^^^iuj, q) + i (^N{lu) + Pl/t(w, q). 



B. The relativistic Langevin equation 

In Sec.[2]we discussed the non-relativistic Langevin dynamics of a heavy-quark in a hot 
plasma. This allowed us to take advantage of well known results like the Einstein relation 
given in Eq. (fTO| , which provides a link between the friction and the momentum-diffusion 
coefficients and can be seen as an example of the fluctuation-dissipation theorem. 

The generalization of the Langevin equation to the relativistic case is far from trivial. 
In particular one can no longer ignore the momentum dependence of the transport co- 
efficients, as it is usually done in non-relativistic studies; in turn, this introduces in the 
definition of the friction term a dependence on the adopted discretization scheme. 

From the mathematical point of view the relativistic Langevin equation belongs to 
the class of stochastic differential equations with multiplicative noise, which represents by 
itself an active field of investigation ^] . Although the traditional Langevin treatment of 
non-equilibrium dynamics is textbook material, the case of state-dependent (momentum 
dependent in the present case ) transp ort coefficients has been only rarely considered in 
the literature [1^, [13, [H, 0, 5^, 51 , and from quite different perspectives. Hence 



we provide in this Appendix a self-contained discussion of the essential aspects of the 
relativistic Langevin equation. We mainly follow the approach presented in Ref. [s^ . 
which we found particularly clear. 

For simplicity, let us consider the one-dimensional case, which is sufficient to illustrate 
the relevant conceptual points. The generalization to the realistic d-dimensional case will 
follow immediately. Also in the relativistic case one postulates that the momentum of 
an external particle placed in a hot medium evolves according to the equation 

^ = -VDip)p + m, (79) 

where the right hand side is given again by the sum of a friction force and a term related 
to the random momentum kicks received from the medium particles 



mm)^4p)s{t-t'), 
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(80) 



still assumed to be uncorrelated. This is expressed by the delta function, leading to a flat 
power spectrum referred to as white noise; however the strength k of the noise in Eq. (I80p 
depends, in the general case, on the momentum of the brownian particle (multiplicative 
noise). One can factor out the momentum dependence by defining g(p) = ^ so 
that Eq. ((75|l becomes 

^ = -'?i3(p)p + gb)'7W, (81) 

with 

m^[t;))^&(t~t;). (82) 

We anticipate that the expression of the friction coefficient rjuij)) will depend on the 
discretization scheme, namely on the value of p at which in the discretized equation 
(employed in the numerical simulations) 

Ap = p{t + At) - p{t) = -[T^D{p)p]At + g{p)v{t)At (83) 

the terms in the right hand side must be evaluated. In the following we wish to clarify 
this point. 

The Langevin equation can be written in the general form 

^^f{p)+9{p)vit), (84) 

where we can identify 

f{p) = -Vd{p)p ^ -Vd^ ip)P + flip) 

p. (85) 



(0), ^ hip) 

Vd ip) 



p 

In the above the friction coefficient has been written as the sum of a leading term, 
linked to the momentum diffusion coefficient k(p) by a relation analogous to the one 
occurring in the non-rclativistic case, and a subleading (in T/E) correction which will 
depend on the discretization procedure. The ambiguity will be eliminated by requiring 
that in the continuum limit all the different discretization schemes yield to the same 
Fokkcr-Planck equation for the momentum distribution, which admits the equilibrium 
distribution exp (^Ep/T) as a steady solution. 

We start by integrating Eq. (|84|). which leads to the formal expression 



/t+At 
ds[fipis))+gipis))Tjis)]. (86) 

However the random noise term 77(5) makes the value of the integral ambiguous and a 
recipe has to be given to evaluate the right hand side of the above equation. One can 
overcome this difficulty by considering a whole family of different discretizations, labeled 
by a parameter a G [0, 1], such that 

l-t+At 

Ap = f[p{t) + aAp]At + g[p{t) + aAp] J dsr]{s). (87) 
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Note that the momentum increment Ap arises from the sum of a deterministic friction 
term of order At and of a random term of order V At. After defining po = the first 
term in the above, being aheady of order Ai, can be evaluated at pq, while for the second 
term it is useful to expand 

g[po + aAp] = g{po) + g'{po) aAp + . . . , (88) 

where the last term, multiplying the noise-integral, provides also a contribution of order 
At. Taking advantage of the fact that {ri{s)) = and of Eq. (f82|) one immediately gets: 



{Ap) = f{po)At + ag{po)g'{po)At 



(89) 



and 



((Ap)2) = g^p,)At. 



(90) 



We now consider the link between the Langevin equation and the Fokker-Planck equation 
for the momentum distribution P{p, t) of the brownian particles. For the latter the 
following equation holds: 



/ + 00 
dpnP{p,t + At\pn,t)P{po,t), 
-OC 



(91) 



where P{p,t + At\pQ,t) represents the conditional probability that a particle with mo- 
mentum pq at time t will be found with momentum p at time t + At. One can identify 
such a conditional probability with the following expectation value over the ensemble of 
brownian particles: 



Pip,t + At\po,t) = {S[p-pit + At)])p„,t 



(92) 



with Ap taken from the Langevin equation (|87p . One can then expand up to second 
order, obtaining: 



P{p, t + At\po,t) = 6{p - po) - (Ap) ^Sip - po) + ^{iApf)^6{p - po) 



(93) 



After inserting the above expansion into Eq. (f9T|) and exploiting Eqs. (l89|) and ((90|) . one 
arrives to the Fokker-Planck equation: 



d_ 

di 



PiP,t) 



d_ 
dp 
d_ 
dp 
d_ 
dp 



-fip)- ag{p)g'{p) 



1 d 



2 dp 

VD'iP)P - flip) ~ o:g{p)g'{p) 

(0) 



1 d 



2 dp 



9\P) 



u 1 Id 

VD'ip)p- hip) + ^ + 



Pip,t) 

P{p,t). (94) 



Then, the requirement that the above equation is independent of the discretization 
scheme (i.e. of a) and admits the steady solution exp i-Ep/T), with Ep = ^p^ -V M2, 
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allows to fix in an unambiguous way rj^jP (p) and fi{p). One gets: 

V^\P) = 1^ (95) 
flip) = ^{l-a)dpK{p). (96) 

Hence the friction coefficient to be employed in tlie numerical Langevin simulations de- 
pends on the discretization scheme through the parameter a and reads: 

Vd{p) ^ - 7;i^ - a)— ■ (97) 

21 hp L p 

Two very popular choices in the literature are a = (Ito discretization) and a — Xjl 
(Stratonovich discretization). For our purposes it will result more convenient to keep 
track of the dependence on the velocity of the brownian particle v = p/ Ep rather then 
on its momentum, so that: 

= (98) 

We notice that the Fokker-Planck equation turns out to be completely determined by 
the momentum-diffusion coefficient: 



p d 



TEp dp 



P{p,t)\. (99) 



We now consider the relativistic brownian motion in d space dimensions. The Langevin 
equation in this case reads: 

^ = -77z,(p)p*-t-f(0, (100) 

with 

{C{t)£,'(t'))^b^^{p)5{t-t'), (101) 

where 

b'^{p) = KL{p)ffP + kt{p){5'^ - fp>). (102) 
It is also useful to introduce the tensor 



g'\p) EE Vkl(p)j5'p^' + ^/H^T{p){5''' -p'p') 

= gL{p)fp' +gT{p){5'' -py). (103) 

This allows one to factor out the momentum dependence of the noise term in Eq. pOOp . 
which becomes 

-^ = -VD{p)f + g'Hp)v'it), (104) 

with 

{r]Ht)r]^{t')) = S'^S{t~t'). (105) 
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Eq. (|104p can be viewed as a particular case of the generic stochastic equation 



(106) 



One can then repeat the same steps foUowed in the one-dimensional case. From the 
Langevin equation one gets the expectation values: 



(Ap') = ripo)At + a (dkg'^po)) ff'^(Po)Ai 
{Ap^Ap^)^g^Hpo)9"'iPo)At = b'^ipo)At, 



which lead to the Fokker-Planck equation: 



d_ 
dt 



P{p.t) 



d_ 
d_ 



1 d 



-f{p) - a {dkg'Hp)) g^'ip) + 2d^b^'^P^ 



P{p,t) 



-fip) + {b^^ip)) - a (dkg'^p)) g'^p) + ^b^p)^ 



(107) 



(108) 



Pip,t). 



Also in the d-dimensional case, requiring the Fokker-Planck equation to be independent 
on the discretization scheme and to admit a relativistic Maxwell distribution as a steady 
solution allows to fix the friction term entering into the Langevin equation and to relate 
it to the momentum-diffusion coefficients. One gets: 



r(p) 

which leads to 

fiP) - 



2T dpi 2 dpi 



{dkg'Hp))9'''ip) 



(109) 



TP 



2TE 
a gL{p)dpgL{p) 



dpKLip) + - — -{kl{p) - kt{p)) 
P 



d-1 
P 



gT{p){9L{p) - grip)) 



(110) 



The friction coefficient to be employed in the Langevin equation can be then conveniently 
written (for the sake of simplicity we give it for the a = Ito discretization) as: 



v'^Tip) 



kl{p) _ J_ 
2TE 



2ydKL{p) , d-l k.l{p) - ht{p) 



(111) 



to be compared with the result given in Rcf. 26]. The recipe employed to update the 
heavy-quark momentum is then the following: 



with 



Pn+l-Pn = -VD°iPn)pl.At + e{tn)At 

= -VD°{Pn)pl,At + g'i{p^)C{tn)VAi, 



(112) 



(113) 



Hence, at each time-step and for each quark, one has simply to extract d independent 
random numbers from a gaussian distribution with cr = 1, as it can be seen from Eq. pi3p . 
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C. Estimate of q, 



In order to obtain a realistic estimate of gmax to employ in the numerical calculations 
we adopt the following strategy. We take, as the upper bound of integration over the 
exchanged momenta, the ultraviolet cutoff given in Ref. [s^. The latter, in the limit 
T required by our eikonal approach to be reliable, reads: 



2(k){E+p) 



^AP + 2{k){E + p) 



(114) 



In the above (fc) ~ T represents an average momentum of the thermal particles taking 
part to the collisions. In Ref. 35| the authors let it vary from T to ST. In order to 
avoid plotting too many curves, our procedure is instead the following. We consider the 
non-relativistic limit (i.e. M 3> T) and we match the results for the momentum diffusion 
coefficient k given by the QCD kinetic calculation [5^ 



3k" 



9 

27r3 



/ 



k^dk 



2k 



dq 



(115) 



2 (e^'= + l)2 
and by the HTL effective approach, namely 



''Dl 

N, 



2 (e/3''— 1)2 



fc2 



4fc4 



3k"'' = Cf 



2tt 



q^dq 



{q^ + ?7i|,)2 ' 



where 



2 2^2 l^c ^Nf 



(116) 



(117) 



We found that, for the case of a coupling running with the temperature, the expression 

Cax - 3.1Tgi/3(T) (118) 

provides a good matching, with an agreement at the 1% level over the whole range of 
temperatures. Inserting the above expression into the zero momentum limit of Eq. (|114p 
allows then to fix {k) , and hence to determine the corresponding bound for general values 
of p. 
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